Orbital mixing and nesting in the bilayer manganites La2-2a:Sri+2xMn207 
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A first principles study of La2-2iSri+2i:Mn207 compounds for doping levels 0.3 < a; < 0.5 
shows that the low energy electronic structure of the majority spin carriers is determined by strong 
momentum dependent interactions between the Mn Cg d^2_y2 and d^,,2_j.2 orbitals, which in addition 
to an X dependent Jahn- Teller distortion, differ in the ferromagnetic and antiferromagnetic phases. 
The Fermi surface exhibits nesting behavior that is reflected by peaks in the static susceptibility, 
whose positions as a function of momentum have a non-trivial dependence on x. 

PACS numbers: 71.18.+y, 71.20.-b, 75.40.Gb, 75.47.Lx 



One of the most studied series of compounds over the 
past decade, having being scrutinized using a wide range 
of experimental methods [l|] , is that of the bilayer man- 
ganites La2-2a:Sri+2a;Mn207, with new findings of their 
basic properties continuing to emerge [3, l3|. Their in- 
terest is due to several reasons. It is well known that 
they exhibit a colossal magnetoresistive effect around 
critical temperatures of 100-200 K, depending on the 
doping level a; [l|. A proper understanding of this phe- 
nomenon should point the way to tailoring materials with 
higher critical temperatures, which could potentially lead 
to magnetic devices that would outperform present ones 
based on the giant magnetoresistive effect. On the fun- 
damental side, it is thought that the properties of these 
materials are due to an interplay of structural, orbital, 
and magnetic degrees of freedom, a trait shared by sev- 
eral strongly correlated materials, including the super- 
conducting cuprates [J| . In the case of the manganites, 
there is no consensus as yet on exactly how these differ- 
ent degrees of freedom lead to the observed properties. 
Elucidating this will no doubt lead to an advance in our 
understanding of related phenomena in the broader con- 
text of higher complexity condensed matter systems. 

One of the reasons for the current state of affairs is that 
the ground state properties of the bilayer manganites are 
not completely understood, even at low temperatures. 
From the theoretical standpoint, there are few published 
reports of ab initio calculations of the electronic struc- 
ture of these compounds [5]. While these studies have 
provided an important framework for the analysis of ex- 
periments, several observations remain to be addressed. 
For instance, it is generally recognized that it is the char- 
acter of the occupied Mn Cg states that determines the 
key properties of these systems, such as the conductivity 
and the nature of the magnetic order. Thus, important 
efforts have been made to understand the orbital polar- 
ization of these states [a lD, [sj, |9|], but the problem has 
only been partially addressed from an ab initio perspec- 
tive [lOj . In this work, we present a first principles study 
of the conducting states of the bilayer manganites for hole 



(i322_r2 orbitals strongly mix, with the mixing depending 
not only on x, but also on the direction in k-space. This 
has a pronounced impact on the Fermi surface topology 
and nesting, long suggested to play an important role in 



these materials [11|, |12|, [13| . This is illustrated by calcu- 
lating the static susceptibility x(q)j which indeed shows 
peaks due to nesting. The implications of our work in 
relation to recent experiments are discussed. 

We employ the highly precise all-electron full-potential 
linearized augmented plane wave (FLAPW) implemen- 
tation of density functional theory [lj|, with the gen- 
eralized gradient approximation (GGA) [15|, ll6| for the 
exchange-correlation potential. For x < 0.5, we use the 
"charged cell" approximation, i.e., increasing the valence 
electron number by the required amount and adding a 
uniform positive charge for neutrality 17|. The calcula- 
tions were carried out for the observed magnetic phases: 
ferromagnetic (FM) for x < 0.4 and antiferromagnetic 
(AF) for X > 0.45. In our calculations, the AF phase 
is given by ferromagnetic Mn02 planes that are coupled 
antiferromagnetically within a bilayer as well as between 
adjacent bilayers. Though in experiment the latter cou- 
pling is ferromagnetic in the A-type AF order, it is known 
to be very weak [l|, and so should not influence our re- 
sults [16|. As for technical aspects of the calculations, 
convergence was assured with respect to muffin-tin radii, 
k-point mesh [1164 points in the irreducible wedge of the 
Brillouin zone], and energy cut-offs. The structural pa- 
rameters for each x (I4/mmm) are from Ref. [l8| at 10 K. 
We assume La is in the 26 and Sr in the 4e sites |19| . 

As indicated above, key to the understanding of the 
bilayer manganites is the relative role of the Cg states 
{dx2-y2 and d^z^-r^) [l|. To illustrate these states, we 
refer to the band structure plots in Fig. [T] (simple tetrag- 
onal symmetry notation is used). Two cases are shown, 
namely x = 0.5 (AF) and x = 0.4 (FM). The states with 
more than 50% Mn d character are highlighted with dots. 
In the AF case, there is Kramers degeneracy, so the two 
spin states are interchanged for the second Mn atom type 
(there is only one Mn atom type in the FM case). In the 



doping levels 0.3 < x < 0.5. We show that the dj.2_y2 and following, we focus on the eg states. 
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FIG. 1: (color online) Band structure of LaSr2Mn207 {x — 
0.5, AF phase) for t [|] spin in panel (a) [(b)], and of 
Lai.2Sri.8Mn207 {x — 0.4, FM phase) for j [J,] spin in panel 
(c) [(d)]. The gray (red) and black (blue) dots highlight states 
of dominant Mn t2g and eg character, respectively. 



The basic electronic structure near the Fermi energy 
(E-p) can be understood as arising from the two Mn 
dj.2_y2 and two Mn d^^2_j,2 orbitals per bilayer. Ignoring 
orbital mixing and hybridization with other states, in the 
FM case the energy of each orbital would be of the form 
Ei ± X ± A, where Ci is the unpolarized atomic orbital 
energy {i = 1,2), 2X the exchange splitting, and 2 A the 
bilayer splitting, resulting in eight states altogether. In 
contrast, in the AF case the energy of each orbital would 
be of the form e^ ± Va^ + A^, resulting in four Kramers 
degenerate states. Since X 3> A, the bilayer splitting is 
essentially quenched in the AF case. Consequently, as we 
show below, there is only one barrel centered at M in the 
AF case, as opposed to two in the FM case. Moreover, 
for the AF case, unlike the FM one, we do not expect a 
strong intensity modulation of the photoemission signal 
as a function of photon energy due to bilayer splitting. 

In the actual calculations the two eg orbitals are mixed. 
This is directly reflected in the Fermi surface topology. 
In Fig.[2l[a), we show the Fermi surface for x = 0.5 (AF), 
and in Fig. ^h) for x = 0.4 (FM, majority spin). In 
the AF case, there is a hole-like barrel centered at M of 
mainly dx2_y2 symmetry (labeled 1), a prolate electron- 
like sheet around F of d^z'^_r2 symmetry (labeled 2), and 
a (minority spin) cylindrical sheet centered at F of dxy 
symmetry (labeled 3). In the FM phase, there is a hole- 
like barrel of dominant dj.2_y2 symmetry centered at M 
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FIG. 2: (color online) Fermi surface for (a) x = 0.5 (AF) and 
(b) X = 0.4 (FM). (c), (d), and (e): Orbital character of the 
states along the F-X-M line for the eg bands crossing Ey for 
X = 0.4. (f) Luttinger count of the different Fermi surface 
sheets, and (g) integrated density of states (IDOS) for d^2_y2 
and d^^2_^2 as a function of doping, imposing either AF or 
FM order. Note the crossover in orbital population due to 
the change in magnetic order (dotted lines). 



symmetry: one hole-like, centered at M (labeled 1'), and 
the other electron-like, centered at F (labeled 2) [20|. 

The mixing is shown by an analysis of the character 
of the states for bands crossing E-p. For example, in 
Fig. m^c) we show the character of the lower tg band (!') 
for X = 0.4 (FM) and k points along the F-X-M direc- 
tions. It starts with dj.2_y2 character at F, but when 
reaching X, the mixing with d^z'^-^'^ is about 50%. At 
kp, roughly halfway between X and M, the character be- 
comes dominantly d^z'^_r2. The middle band (1), in con- 
trast [Fig. [2)Jd)], has a fairly dominant dr^2_y2 character 
up to fcp, at about a quarter of the X-M distance. Band 
2, in turn [Fig. [D^e)], starts out at F as d^z^_j.2 but is 
strongly mixed at fcp. Thus, it is inaccurate to charac- 
terize the barrels around the M point as 'x^ — y^\ and 
the electron pocket around F as '3z^ — r^\ Our results 
for the other doping levels in the FM phase show that 
the orbital mixing increases with decreasing x, and that 
the bilayer splitting becomes stronger. A similar analysis 
shows that in the AF case, the lower Cg band (1) is of 



dominant d,^ 



character, but with strong mixing at kp 



(labeled 1), and two barrels of mixed d^2_y2 and d; 



3z2-r-' 



along X-M. The upper band (2) is d^z'^_j.2 like. 

Of further relevance to experiment is the occupation 
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FIG. 3: (color online) The static susceptibility, x(ci)j for Q ^ 
{q^, qy,2n/c), for x = 0.5 (AF) and x = 0.4 (FM). 



level of the bands crossing Ep. Thus, we calculated the 
carrier number by integrating the volume enclosed by the 
different Fermi surfaces (the Luttinger count). We used 
a fine mesh of 32340 k-points in the irreducible wedge, 
and the results as a function of x are shown in Fig.[2lf). 
In the AF phase, sheet 1 is obviously dominant. In the 
FM phase, the Kramers degeneracy is lifted, giving rise 
to sheets 1 and 1'. The increased orbital mixing lowers 
the energy, particularly in the case of states 1', allow- 
ing the charge to increase as x decreases. In the FM 
phase there is only one sheet 2 because the antibonding 
counterpart is empty (due to the large bilayer splitting 
of 3z^ — r^ in the FM phase). Thus, the electron pocket 
around F roughly doubles in size compared to the AF 
phase. Note that for 0.4 < x ^ 0.45 the actual phase 
is a canted antiferromagnet [18|, so the bilayer splitting 
is probably reduced gradually with x. But clearly, AF 
coupling acts to quench the bilayer splitting, with dra- 
matic effects, particularly for the d^z2_j,2 states. This is 
in line with very recent photoemission data that indicate 
a collapse of bilayer splitting near a;=0.4 2lU22|. 

In Fig. [2Ig) we plot the Cg orbital occupation counts 



(estimated by integrating the muffin-tin projected den- 
sity of states) as a function of x, imposing either a FM or 
AF phase. The d^^^y^ orbital population is not strongly 
affected across the FM-AF transition, but the d3z'^_r^ or- 
bital population significantly increases in the FM phase. 
This agrees with conclusions from magnetic Compton 
scattering measurements [6|. Further, our results show 
a striking cooperative effect between Jahn- Teller distor- 
tion and magnetic order. Indeed, the evolution of the oc- 
cupations with X in a given magnetic phase reflects the 
change of the apical Mn-0 bond lengths [l^l- But the 
crossover due to the change in magnetic order (dotted 
line) is largely due to the change in bilayer splitting. 

The above results also indicate that the possible nest- 
ing instabilities s ugg ested to play an important role in 
these materials [ll|, ll3| have a non-trivial dependence 
on x. To see this, we calculated the generalized charge 
susceptibility (constant matrix element approximation) 
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FIG. 4: (color online) Static susceptibility x(ci) for Q along 
two symmetry directions for x — 0.5 (AF) and x — 0.4 (FM). 
The labels indicate which band transitions are responsible for 
the structure in x- Labeling corresponds to that in Fig. [S] 
Note the break in the vertical axis. 



For this calculation, we again used 32340 k-points in the 
irreducible wedge and the tetrahedron method with lin- 
ear interpolation [23|. In Fig. [31 we show x(q) for q of 
the form (q^, qy, 2tt/c). This is of interest because diffuse 
x-ray scattering data have shown peaks for q^ = 2tt/c 
in the x 



0.4 compound [2j|. Both the AF and FM 
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cases show clear peaks and kinks correlating closely with 
the topology of the Fermi surface. The calculations were 
done taking into account only the bands crossing Ep . We 
verified that including more bands does not affect signif- 
icantly the momentum structure of the response. Also, 
the Cg to t2g transitions were ignored to crudely simu- 
late the neglected matrix elements effects. We note that 
the t2g transitions mostly contribute to a diffuse response 
centered at T, so we focus on the eg transitions only. 

Figure [4] takes a closer look for q along two directions. 
For q = {q, 0, 2tt/c), Fig. \^a,) (AF) shows a strong peak 
at g = 0.53 (11 transition). The interband transition 
(12) has little structure (not shown). Figure IDJb) (FM) 
shows two main peaks, for q = 0.62 (I'l transition), and 
0.71 (I'l' and 22 transitions). The peak arising from 
11 transitions is relatively weaker in the FM case. It is 
significant that the I'l peak at this doping level matches 
the bond centered peak reported at q = 0.6 in Ref. ]2d- 
For q = {q,q,2'K/c), Fig. gj^c) (AF) shows the same (11) 
peak as Fig. Ufa), and a less well defined maximum at 
q = 0.42 (12 transition). This maximum is equally ill- 
defined at other doping levels. Fig. \MA) (FM) shows a 
clearly dominant peak at g = 0.7 arising from mainly I'l' 
and 1'2 transitions. We point out that the two kinks for 
smaller q (I'l and 11) are weaker at lower doping levels 
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FIG. 5: (color online) Nesting vector lengths (extracted from 
the peaks in x) ^-s a function of doping level x for q along (a) 
the [100] and (b) [110] directions. 



the Cg orbital polarization, are predicted to have a non- 
monotonic dependence on doping due to the change from 
FM to AF order as the hole doping level increases. 
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because of strong broadening of these features. 

Finally, in Fig. [5] we consider the dependence on x 
of the lengths of the magnitude of the nesting vectors 
along the [100] and [110 ] directions, inferred previously 
from photocmission [ll|, [l^. We extracted these from 
the peaks in x corresponding to those in Fig. |31 Again, 
we remark that we do not have results for the actual 
canted phase for 0.45 ^ x < 0.5. Note that the 11 vec- 
tor for X — 0.4 in the [110] direction matches that of 
the so-called CE-ordering wave vector jll|. We partic- 
ularly point to the jump of the 11 nesting vector when 
the magnetic phase changes. Clearly, this is due to the 
'turning on' of the bilayer splitting in the FM phase. 
Again, this illustrates dramatically the correlation be- 
tween d^z^-r^ and dx^-y^ orbital mixing, bilayer split- 
ting, and FM order. Further, this jump may explain 
the non-monotonic x-dependence of the nesting vector re- 
cently inferred from photoemission and diffuse scattering 
measurements [25|. The charge ordering peaks observed 
in the latter are either along the bond direction, the diag- 
onal, or both, depending on the doping level. Typically, 
the diagonal response is largest, since although the nest- 
ing is not as good for such vectors, twice as many Fermi 
surface faces are brought into coincidence as compared 
to the bond oriented case [26[. We also note that in the 
FM phase, for transitions involving 1', the peak structure 
for Qz = along the bond direction is less distinct than 
for Qz — l-KJc. This is because of the 'tilting' of the 1' 
surface (a consequence of its strong d-^z'^-^'^ admixture), 
which alternates in direction between gz = and 27r/c 
due to the c-axis dispersion. For the same reason, the 
diagonal peaks involving 1' are stronger for qz = 0. This 
may explain why the diffuse scattering peaks occur for 
qz = 27r/c for the bond directions and ^z = for the 
diagonal directions 2J, |25| . 



In summary, we find that mixing of the Mn (i^2_j,2 and 
d322_^2 orbitals and bilayer splitting play a fundamental 
role in the electronic structure of the bilayer manganites, 
particularly in the FM phase, and correlate closely with 
the transition from FM to AF order with doping. These 
two effects result in different Fermi surface topologies for 
the two phases. As a result, the static susceptibility, and 
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